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Abstract In this paper we review the present status and impUcations of Ha sur- 
veys at various redshifts. With the advent of sensitive wide-format 
near-infrared detectors on large telescopes, deep and extensive Ha sur- 
veys are now feasible to redshift z ~ 2.5. The cosmic star formation 
history can therefore be traced out to this redshift using Ha alone, 
avoiding complications arising from the comparison of different tracers 
at different redshifts. The Ha surveys to date confirm the rapid in- 
crease in luminosity density from z — out to 2: — 1, and show that 
this increase flattens off at higher redshifts, remaining approximately 
constant out to at least z ~ 2.2. We also discuss the prospects for de- 
termining the masses of high redshift galaxies based on emission lines. 
A set of high-quality Ha rotation curves of samples of disk galaxies at 
a number of different redshifts would allow a study of the evolution of 
the TuUy-Fisher relation and address fundamental issues in disk galaxy 
formation. Such a program remains challenging even with present-day 
large telescopes. 
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1. Introduction 

When did the stehar populations of present-day galaxies form? This 
question is now beginning to be addressed by observations with mod- 
ern telescopes in virtually all regions of the electromagnetic spectrum. 
Moderate redshift (z < 1) surveys consistently indicate a cosmic star 
formation rate density (SFRD, the stellar mass formed per unit of time 
and per unit of comoving volume) that strongly increases with red- 
shift out to z ~ 1 (Hogg, 2001). The SFRD at z > 1 is however a 
more hotly debated issue. Estimates based on surveys of Lyman break 
galaxies (e.g., Madau et al., 1996) sample rest- frame UV radiation and 
require a large and uncertain extinction correction (e.g., Pettini et al., 
1998; Steidel et al., 1999). Submillimetre surveys, tracing reradiation 
by dust of absorbed starlight (e.g., Hughes et al., 1998; Blain et al., 
1999) consistently indicate an even higher SFRD at these redshifts than 
even the extinction-corrected UV-determined values. As a result, the 
use of different star formation tracers at different redshifts complicates 
the reconstruction of the cosmic star formation history and it is there- 
fore highly advantageous to use the same star formation tracer at all 
redshifts. The Ha line is the natural choice for this tracer, since with 
the development of large area, high quality near-infrared arrays, this 
line can now be observed in the J, H and if-bands, where redshifts of 
approximately 1, 1.5 and 2.2 respectively, are accessible. 

Traditionally, emission-line surveys for high-2; star forming galaxies 
have targeted the Lya line, which moves into the optical regime for 
z = 1.9 — 7.0. However, since Lya is resonantly scattered, even very 
small quantities of dust will effectively suppress the line (Chariot & Fall 
1991, 1993; Chen and Neufeld, 1994). Spectroscopic observations of star 
forming z > 3 Lyman break galaxies (Steidel et al., 1996) confirm that 
Lya is not a good tracer of star formation in high-z galaxies. These 
galaxies form stars at rates of ~ 10 Mq yr~^ as derived from their 
restframe UV properties (uncorrected for extinction), but Lya is absent 
(or in absorption) in more than 50% of the cases, while in most of the 
remaining objects the line is faint. Thus, while Lya searches can be 
used to find high- 2: star forming galaxies (e.g., Hu and McMahon, 1996; 
Hu et al., 1998; Cowie and Hu, 1998; Manning et al., 2000; Stern et al., 
2000; Kudritzki et al., 2000; Rhoads et al., 2000; Steidel et al., 2000), 
such surveys cannot provide reliable star formation rates (SFRs) for the 
galaxies that are detected. 

This complication is avoided by searching for Ha emission in stead 
of Lya. The Ha line is not resonantly scattered and thus much less 
sensitive to the effects of small amounts of dust. In addition, the broad- 
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band extinction at Ha (6563 A) is much less than that at Lya (1215 A) 
(^LyaMna = 4.28 for the extinction curve used by Cardelh et al., 1989). 
Although the extinction at the wavelength of Ha is still large in star- 
burst galaxies, the typical Ha extinction in local spiral galaxies is only 
1™1 (Kennicutt, 1983 and references therein). A recent comparison of 
star formation tracers (including stellar continuum emission as well as 
nebular lines) confirms that Ha is the most robust star formation tracer, 
and can give results accurate to within a factor of 2 to 3 when combined 
with other tracers, even if the Ha/ [Nil] blend is not resolved (Chariot 
and Longhetti, 2001). In contrast, the also widely used [Oil] line at 
3727 A is a more uncertain star formation tracer by a full order of mag- 
nitude, because of uncertainties in excitation, metallicity and reddening. 
The SFR derived from Ha also correlates well with that estimated from 
the UV (Bell and Kennicutt, 2001). The systematically higher SFR de- 
rived from Ha (by a factor of 1.5 in relatively luminous galaxies) simply 
reflects the difference in extinction between the optical and UV regimes. 

In this paper we first review the current status of Ha surveys from 
z = out to the highest redshifts where such surveys have been per- 
formed to date. Since we only consider Ha, this analysis should be free 
of uncertainties introduced by using different tracers at different red- 
shifts. Uncertainties in the conversion of Ha luminosity into SFR rate 
can be avoided by quoting the Ha luminosity density, which can be de- 
rived directly by integrating the Ha luminosity function. 

In the second part of this review we address the use of Ha for study- 
ing the kinematics of high-z galaxies. While the SFRD is expressed as a 
time derivative of comoving mass density, the various ways of determin- 
ing SFRD do not involve the measurement of mass but of light. However, 
the stars producing most of the light contain very little mass, while the 
stars containing most of the mass produce little light, for a typical ini- 
tial mass function (IMF). As a result, calculating the total SFR from a 
measurement of light (regardless of whether this is rest-frame ultravio- 
let light, nebular line emission, or emission from dust) always involves 
the uncertain step of assuming an IMF. While systematic errors in the 
analysis can be mitigated by ensuring that always the same IMF is as- 
sumed, the universality of the IMF is still an assumption and a very 
debatable issue. The most direct way of avoiding this complication is 
by making a dynamical mass measurement. Spectral line measurements 
offer this possibility, provided the kinematic pattern is dominated by 
gravitational motion. For this reason, Lya kinematics is normally not 
useful since Lya line profiles are often affected by nongravitational mo- 
tions such as outflows. Emission line kinematics can be assumed to be 
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dominated by gravitational motion if a characteristic pattern such as a 
rotation curve is observed. Even then, one must be wary of excitation 
effects, since Hnes such as [O III] 5007 A possibly sample only the inner 
parts of a galaxy, and may not reveal the full amplitude of the rotation 
curve. Again, using Ha is the best way for avoiding such biases. 

Dynamical measurements of this type have the potential for probing 
the evolution of the Tully-Fisher (TF) luminosity-line width relation 
out to z ~ 2.5. This issue is directly related to the physics of disk 
galaxy formation. One model for the origin of the TF relation argues 
that the dark halo is the sole factor determining the disk properties 
(and therefore the TF relation), since the outer region of disks, where 
the rotation curves are flat, are dynamically dominated by the dark 
halo. Models based on this observation predict rapid evolution of the 
TF relation (Mo et al., 1998). On the other hand, it has been proposed 
that the TF relation is entirely the result of star formation self-regulation 
in the disk interstellar medium and has nothing to do with the halo at 
all (Silk, 1997). A set of high-quality Ha rotation curves of disk galaxies 
at a number of different redshifts is needed to settle this fundamental 
issue from an observational point of view. First results and prospects in 
this challenging field are reviewed in Sect. 3. 

Unless noted otherwise, in this review, all quantities that are de- 
pendent on cosmology are quoted for a universe with {h,^lM,^A) = 
(0.5,0.5,0.0), where the Hubble constant is Hq = 100 h km Mpc"^ 
and JIm and i}\ are the present-day density parameters of matter and of 
the cosmological constant, respectively. Where necessary, results origi- 
nally quoted for other cosmologies have been tacitly converted. 

2. The cosmic star formation history derived 
from HcK surveys 

2.1. Surveys and strategies 

Observationally, the task of measuring the SFRD at a particular red- 
shift from Ha consists of determining the Ha luminosity function, which 
is parametrized (Schechter, 1976) as 

(t){L) dL = (j)* {L/L*)'^ e-^/^' d{L/L*) (1) 

and gives the number density (per Mpc^) of galaxies with luminosities 
between L and L + dL. For cosmological applications, comoving number 
densities (i.e., number of galaxies per comoving Mpc^) are computed. 
In practice, the number density of galaxies is binned per 0.4 interval in 
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log L, which can be related to the standard form by 




(2) 



The total luminosity density is found by integrating the luminosity func- 
tion over all luminosities: 



where T denotes the gamma function. For a determination of the to- 
tal luminosity density it is therefore sufficient to determine the three 
parameters (j)*, L* and a of the luminosity function. Note that the in- 
tegrated luminosity density is dominated by low- luminosity objects if a 
approaches —2. 

Surveys for emission line objects can be of three types (cf., Djorgov- 
ski, 1992). Slitlcss spectroscopic surveys sample large volumes of space 
but can only achieve interesting depths when carried out from space. 
Ground-based surveys can be either serendipitous long-slit surveys (cov- 
ering a large redshift range but a very limited solid angle) or narrow-band 
surveys (sampling a large area but in a very narrow redshift interval). 
Since the background in the near-IR windows shortward of about 2.2 /zm 
is dominated by OH emission lines, the narrow-band STirvcy technique 
is the ideal approach in this spectral region for ground-based observa- 
tions. This technique involves deep imaging in a suitable narrow-band 
filter, complemented with broad-band imaging; sources with excess flux 
in the narrow-band filter are emission line candidates in a redshift in- 
terval determined by the narrow-band filter passband (this strategy has 
been analyzed in detail by Mannucci and Beckwith, 1995). Narrow-band 
surveys of Ha emission at 2 > 2 were initially unsuccesful (Thompson 
et al., 1994; Pahre and Djorgovski, 1995; Bunker ct al., 1995; Collins 
et al., 1996; Thompson et al., 1996) due to lack of depth and/or cover- 
age (with the exception of one object at z = 2.43 reported by Beckwith 
et al., 1998). In contrast, more extensive Ha surveys targeted at vol- 
umes containing known damped Lya and mctal-linc absorbers, quasars 
and radio galaxies resulted in the identification of a significant number 
of candidate objects (Mannucci et al., 1998; Teplitz et al., 1998; van der 
Werf et al., 2000). However, since none of these surveys has been suc- 
cesfully followed up with spectroscopic confirmation, these results must 
so far be considered tentative. More fundamentally, since these were 
targeted surveys, centred at particular "marker" objects of known red- 
shift, they are most likely biased to overdense regions and cannot be 
used to derive a star formation rate density that is valid for the universe 
on global scales. 
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Figure la. Narrow-band image at 2.09 /im centred on the HDFS WFPC2 field. The 
field size is 5' x 5' with North at the top and East to the left. The squares identify 
spectroscopically confirmed emission-line sources (from Moorwood et al., 2000). 



2.2. Blank field Ha surveys 

2.2.1 Surveys at z ^ 1. The Ha luminosity density of the 

local [z < 0.045) universe has been determined from the Universidad 
Complutense de Madrid objective prism survey (Gallego et al., 1995). 
At somewhat higher redshift, the Ha luminosity density derived from a 
(z) ~ 0.2 sample taken from the /-band selected Canada France Red- 
shift Survey (CFRS) is about a factor of 2 higher (Tresse and Maddox, 
1998). A small survey of 13 CFRS galaxies at z « 0.9 produced an Ha 
luminosity density approximately 10 times higher than that at z ~ 
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Figure lb. Broad-band Ks image corresponding to the 2.09 nm narrow-band image 
of the HDFS WFPC2 field (from Moorwood et al, 2000). 



(Glazcbrook et al., 1999). All of these results include a small (approxi- 
mately 1 magnitude) correction for extinction, and a correction for the 
effects of the [N ll] lines on the derived Ha flux. The magnitude of this 
increase is in excellent agreement with the increase in cosmic luminosity 
density measured in the rest-frame ultraviolet (Lilly et al., 1996) and in 
radio continuum emission (Haarsma et al., 2000). The redshift region 
out to z ~ 1 is thus characterized by a strong increase in luminosity 
density with lookback-time. While the precise values of the correspond- 
ing star formation rate densities are uncertain due to the complications 
discussed in Sect. 1, the increase out to z ~ 1, when measured in one par- 



8 



ticular tracer, is now well documented and is quantitatively consistent 
between various different tracers (Hogg, 2001). 

It should be emphasized that existing Ha surveys at z ^ 1 are limited 
in both width (leading to poor statistics, especially at the high lumi- 
nosity end) and depth (leading to a poor determination of the faint-end 
slope of the luminosity function). Indeed, a recent deeper narrow-band 
survey for Ha at redshifts of approximately 0.08, 0.2 and 0.4 (Jones and 
Bland-Hawthorn, 2001) suggests the need for an upward revision of the 
local Ha luminosity density, mainly as a result of a higher detection 
rate of faint galaxies. This result is in agreement with an analysis of 
the KPNO International Spectroscopic Survey (Gronwall, 1999), based 
on 1126 Ha emitters at 2; < 0.085 detected using the objective prism 
technique. 

2.2.2 Surveys at 1 ^ 2 ^ 2. The redshift range between 
redshifts 1 and 2 is of fundamental interest since it is in this range that 

the cosmic star formation history starts to turn over from a fairly high 
level of star formation rate density (perhaps to first order constant with 
redshift at z ;^ 1) into the rapid decrease discussed in the previous 
section. The shape, origin, and redshift of this turnover are not well 
known. This lack of knowledge is partly explained by the lack of suitable 
redshift tracers in the optical region for redshifts between 1 and 2. In 
addition, at z ;^ 1 the Ha line has moved into the near-infrared regime 
and sensitive large-format near-infrared cameras and spectrographs at 
large telescopes have only recently become available. 

The most important Ha surveys at 1 ^ 2; ^ 2 to date are the two 
slitless grism surveys carried out with NICMOS on the Hubble Space 
Telescope. The first of these surveys (McCarthy et al., 1999) covered 
approximately 64n' and yielded an Ha luminosity function (Yan et al., 
1999) with L*(Ha) exceeding the z = value (Gallego et al., 1995) by a 
factor of 7, implying significant luminosity evolution. The second survey 
(Hopkins et al., 2000), covered only 4.4n' but extended the Ha lumi- 
nosity function towards fluxes fainter by a factor of 2, pinning down the 
faint-end slope of the luminosity function at a = —1.6. In total these 
surveys produced 70 Ha emitters between z = 0.7 and 1.9. A significant 
fraction of these have been confirmed as Ha emitters by spectroscopy of 
other emission lines. Since the two grism surveys define the Ha lumi- 
nosity function at {z) ~ 1.3 over a significant range in luminosity, the 
implied Ha luminosity density is well-determined. The resulting value is 
3.6 • 10^ L0 Mpc~^, without any correction for extinction. It is remark- 
able that this value is 15% smaller than the z = 0.9 value (Glazebrook 
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et al., 1999), and the steep rise in luminosity density from z = to 1 
therefore appears to flatten off at z > 1. 

2.2.3 Surveys at z > 2. The highest redshift survey for Ha 
emitters used the SOFI near-infrared camera on the ESO New Tech- 
nology Telescope to target a lOOD' area in the region of the Hubble 
Deep Field South (HDFS), in a narrow {Az « 0.04) redshift inter- 
val at z ^ 2.2, using the narrow-band imaging technique (Moorwood 
et al., 2000). This redshift interval was chosen because it redshifts Ha 
to 2.09 /7,m, which is a spectral region relatively free of bright OH lines 
and not strongly affected by the thermal background that becomes dom- 
inant at somewhat longer wavelengths. The importance of this survey 
lies in the fact that it is the first succesful blank-field survey of this type 
and at this redshift, and the only one with substantial spectroscopic 
confirmation. One field from this survey is shown in Fig. 1. It reveals 5 
emission line objects, all of which are spectroscopically confirmed. A set 
of confirmation spectra, obtained with ISAAC at the ESO Very Large 
Telescope, is shown in Fig. 2. The number of detections from this survey 
matches precisely the number expected under the assumption that the 
Ha luminosity function at z ~ 2.2 is identical to that at z ^ 1.3. The 
number density of of Ha emitters at z ~ 2.2 is also comparable to that of 
Lyman break galaxies with similar SFRs at z = 3.0 — 3.5 (Steidel et al., 
1996). These results provide further evidence for a constant SFRD at 
z>l. 

2.3. The evolution of the Ha luminosity density 
and star formation in galaxies 

Comparison of the luminosity functions at 2; ^ (Gallego et al., 1995) 

and z 1.3 (Yan et al., 1999) and higher shows that both luminosity 
evolution and density evolution are implied. However, luminosity evolu- 
tion dominates (a factor 7 increase in L* from 2; = to z = 1.3 compared 
to a factor 2.6 increase in (/)*). Indeed, in the z ~ 2.2 sample the implied 
SFRs are 20 to 35 Mq yr-i (for the Ha/SFR conversion factor of Ken- 
nicutt, 1998, which is appropriate for continuous star formation and a 
Salpeter IMF, at solar metallicity) . These values are significantly higher 
than in typical nearby galaxies, but lower than in extreme starburst 
galaxies, such as are found in submillimetre surveys (e.g., Ivison et al., 
2000). This result suggests that, while the submillimetre measurements 
select extreme starburst objects which may plausibly be identified with 
the formation of present-day spheroids, the Ha measurements select a 
less extreme mode of star formation. It may then be speculated that the 
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Figure 2. ISAAC spectra of emission line galaxies from the HDFS narrow-band 
imaging survey (Moorwood et al., 2000). The markers under the redshift labels show 
the expected wavelengths of the [N ii] (6548 and 6584 A) lines assuming the detected 
line is Ha. 



Ha measurements trace the buildup of disks. Dynamical measurement 
(cf., Sect. 3) may shed light on this hypothesis. 
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Figure 3. The evolution of the cosmic Ha luminosity density (lefthand panel) and 
SFRD (righthand panel, calculated using the Ha/SFR conversion factor of Kennicutt, 
1998). In the lefthand panel, open symbols denote measurements that include an 
extinction correction based on the Balmcr decrement. In the right-hand panel, filled 
symbols denote the Ha-based measurements while open symbols arc based on rest- 
frame UV radiation; the arrows on the 2 high-redshift points in this panel indicate 
the magnitude of the estimated extinction correction. In the righthand panel, the 
righthand vertical axis denotes metal production rate density, which is less sensitve 
to the assumed IMF than star formation rate density (Madau et al., 1996). 



The evolution of the Ha himinosity density and the impUed SFRD 
history are summarized in Fig. 3. The Ha results give a consistent 
picture of evolution out to z ~ 2.2. Without extinction corrections, the 
Ha surveys imply higher SFRDs than the UV-based surveys, but the 
discrepancy is typically only a factor of 2 (Fig. 3). It is interesting that 
the extinction corrections proposed for Lyman break galaxies at higher 
redshifts lead to an SFRD very similar to that implied by the Ha surveys 
(see Fig. 3); it should be noted however, that these corrections are not 
very certain. 

3. Dynamical mass measurements at high 
redshifts using Ha 

The low-redshift TF relation derived from Ha rotation curves is well- 
defined (Courteau, 1997). Attempts to detect evolution in the zero-point 
of the TF relation using Ha measurements of higher redshift samples are 
so far inconclusive: while Ha rotation curves out to z ~ 1 implied only 
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Figure 4- Long-slit Ha spectrum of a galaxy at 2 ~ 2.2 showing the signature of 
ordered rotation (left panel); the implied rotation curve is shown in the right-hand 
panel (Moorwood et al., 2000). 

about 0™4 evolution in B (Vogt et al., 1996, 1997), a study of [Oil] 
rotation curves at (z) = 0.35 suggested luminosity evolution by 1.5 — 2 
magnitudes in B (Simard and Pritchet, 1998). These results may be 
reconciled if in the [O ll] data the rotation curves arc not sampled all the 
way out to the flat part. The discrepancy underlines the need for deep 
uniform Ha rotation curves for projects of this type. 

Since the effects of evolution should be largest at high redshift, Ha 
dynamics is particularly important for the most distant galaxies from 
the surveys described in Sect. 2. One of the galaxies from the z ~ 2.2 
SOFT/ISAAC survey (Moorwood et al., 2000) actually shows a tilt in 
the 2-dimcnsional spectrum (taken along a position angle within 10° of 
the major axis of the galaxy) indicating ordered rotation (Fig. 4). With 
the short integration time of 1 hour only than the central portion of 
the rotation curve is seen, so that only a lower limit to the total line 
width can be determined. The resulting peak-to-peak velocity width is 
247 lb 30 km s"-*^, which increases to 275 it 30 km s"-*^ after applying a 



Probing the evolution of galaxies using redshifted Ha emission 



13 



correction for inclination. Thus a rotational velociy of 138 it 15 km s~ is 
implied at a radius of 3 kpc. These parameters are comparable to those 
of local disk galaxies which have rotation curves flattening at velocities 
between 100 and 300 km s~^ at radii of 1 to 5 kpc. It is therefore evident 
that well-developed, relatively massive systems already occur at 2; ~ 2.2. 
More interestingly, with a rest-frame absolute B magnitude Mb = —22.4 
this object lies about 3 magnitudes above the local B-hand TF-relation. 
This conclusion can only be avoided if the full rotation curves extends to 
over 1000 km s~^ and only flattens at a radius larger than 12 kpc which 
would imply an extremely massive galaxy and is therefore unlikely. 

4. Conclusions and outlook 

Surveys of Ho emission are now beginning to probe the cosmic star 
formation history out to z ~ 2.2 in a meaningful manner. However, 
this field is still in a very early stage of development. A fundamental 
problem is that the samples are small, leading to uncertainty in the de- 
rived luminosity function. This is best illustrated by the fact that even 
at z ^ 0, where the available samples are largest, significant discrep- 
ancies remain (cf., Sect. 2.2.1). Great progress in determining the Ha 
luminosity density at low redshifts should be made as a result of the 
Sloan Digital Sky Survey and associated follow-up. At higher redshifts, 
the VIRMOS survey at the ESO Very Large Telescope should provide 
a census out to redshifts beyond 1. Redshifts z > 2 bring Ha into the 
K-hand and an efficient survey would require a cryogenic multi-object 
spectrograph. At z > 2.5 Ha shifts into the thermal infrared and the 
Next Generation Space Telescope will provide the first opportunity for 
deep Ha measurements at these redshifts. 

As important as the composition of good samples is a good estimate 
of the extinction. Ha extinctions can be determined from the Balmer 
decrement, and should be feasible for most galaxies where Ha can be 
detected. Furthermore, for a reliable determination of star formation 
based on Ha, other spectral features should also be observed (Chariot 
and Longhetti, 2001). 

First results on Ha kinematics at high redshift demonstrate the feasi- 
bility of determining rotation curves of high- 2; galaxies with present-day 
instrumentation. Compiling a sample of high-quality Ha rotation curves 
at a number of redshifts out to z 2.2 will be extremely valuable for our 
understanding of the disk galaxy formation and the origin of the Tully- 
Fisher relation. Such a project should ideally make use of the multiplex- 
ing capabilities of future cryogenic multi-integral-field-unit near-infrared 
spectrographs. However, a first sample of well-measured rotation curves 
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out to z ^ 2.2 can already be built up with deep near-infrared Ha spec- 
troscopy of galaxies from the existing samples. This is within reach with 
present-day near-infrared spectrographs on large telescopes. 
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